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ABSTRACT 

GRB980519 is characterized by its rapidly declining optical and X-ray afterglows. 
Explanations of this behavior include models invoking a dense medium environment 
which makes the shock wave evolve quickly into the sub-relativistic phase, a jet-like 
outflow, and a wind-shaped circumburst medium environment. Recently, Frail et al. 
(1999a) found that the latter two cases are consistent with the radio afterglow of 
this burst. Here, by considering the trans-relativistic shock hydrodynamics, we show 
that the dense medium model can also account for the radio light curve quite well. 
The potential virtue of the dense medium model for GRB980519 is that it implies 
a smaller angular size of the afterglow, which is essential for interpreting the strong 
modulation of the radio light curve. Optical extinction due to the dense medium is 
not important if the prompt optical-UV flash accompanying the 7-ray emission can 
destroy dust by sublimation out to an appreciable distance. Comparisons with some 
other radio afterglows are also discussed. 

Key words: gamma-rays: bursts - hydrodynamics: shock waves 



X 



In the standard model of gamma-ray bursts (GRBs) (see 
Piran 1999 for a review), an afterglow is generally be- 
lieved to be produced by the synchrotron radiation or in- 
verse Compton scattering of the shock-accelerated electrons 
in an ultra-relativistic shock wave expanding in a homoge- 
neous medium. As more and more ambient matter is swept 
up, the shock gradually decelerates while the emission from 
such a shock fades down, dominating at the beginning in X- 
rays and progressively at the optical to radio energy bands 
(Meszaros & Rees 1997; Waxman 1997a; Wijers & Galama 
1999). In general, the light curves of X-ray and optical 
afterglows are expected to exhibit power-law decays (i.e. 
F v oc t~ a ) with the temporal index a in the range 1.1 — 1.4, 
given the energy spectral index of electrons p ~ 2 — 3. The 
observations of the earliest afterglows are in good agreement 
with this simple model (e.g. Wijers, Rees & Meszaros 1997; 
Waxman 1997b) . However, over the past year, we have come 
to recognize a class of GRBs whose afterglows showed light 
curve breaks ( e.g. GRB 990123, GRB 990510; Kulkarni et 



al. 1999a; Harrison et al. 1999) or steeper temporal decays 
(i.e. F u oc t~ 2 ; e.g. GRB 980519, GRB 980326; Bloom et al. 
1999). Explanations for these behaviors include three scenar- 
ios: 1) a jet-like relativistic shock has undergone the tran- 
sition from the spherical-like phase to a sideways-expansion 
phase (Rhoads 1999), as suggested by some authors (e.g. 
Sari, Piran & Halpern 1999; Kulkarni et al. 1999a; Harrison 
et al. 1999). 2) the shock wave propagates in a wind-shaped 
circumburst environment with the number density n oc r~ 2 
( Dai & Lu 1998; Meszaros, Rees & Wijers 1998; Chevalier 
& Li 1999; Chevalier & Li 2000; Li & Chevalier 1999 ); 3) a 
dense medium environment in ~ 10 5 — 10 6 cm -3 ) makes the 
shock wave evolve into the sub-relativistic phase after a short 
relativistic one (Dai & Lu 1999a, b) . In the last model, since 
an afterglow from the shock at the sub-relativistic stage de- 
cays more rapidly than at the relativistic one, we will expect 
a light curve break or a long-term steeper decay, depending 
on the time when it begins to enter into the sub-relativistic 
stage. This scenario has reasonably interpreted the break 
in the R-band afterglow of GRB 990123 (Dai & Lu 1999a) 
and the steep decays of the X-ray and optical afterglows of 
GRB 980519 (Dai & Lu 1999b). 
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Recently, Frail et al. (1999a) tried to test the first two 
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models (the jet and wind cases) by means of the radio af- 
terglow behavior of GRB 980519 and found that the wind 
model described it rather well. Due to the strong modu- 
lation of the light curve, however, they could not draw a 
decisive conclusion for the jet case. In this paper, we will 
examine the possibility of describing the evolution of the ra- 
dio afterglow in terms of the dense medium model. Since this 
scenario involves the transition phase of the shock wave from 
the relativistic stage to the sub-relativistic, we have consid- 
ered the trans-relativistic shock hydrodynamics in the nu- 
merical study. We first present the asymptotic result of the 
fitting of the radio data in section 2, and then the numerical 
result in section 3. In section 4, we show that the optical 
extinction due to the dense circumburst medium is not im- 
portant, since the prompt optical- UV radiation, caused by 
the reverse shock, can destroy the dust by sublimation out 
to a substantial distance, as proposed by Waxman & Drainc 
(1999). Finally, we give our discussions and conclusions. 



2 ASYMPTOTIC BEHAVIOR OF THE RADIO 
AFTERGLOW IN THE SUB-RELATIVISTIC 
STAGE 

GRB 980519 is the second brightest GRB in the Bep- 
poSAX sample. Its optical afterglow measured since ~ 8.5 
hours after the burst exhibited rapid fading, consistent with 
t -2.os±o.04 in BVRI (Halpern et al. 1999; Djorgovski et al. 
1998), and the power-law decay slope of the X-ray afterglow, 
ax = 2.07±0.11 (Owens et al. 1998), was in agreement with 
the optical. The spectrum in optical band alone is well fitted 
by a power-law I/ - 1 - 20 ±°- 28 j while the optical to X-ray spec- 
trum can also be fitted by a single power-law ;/ - 105 ± 01 " 
The radio emission was observed with the Very Large Ar- 
ray (VLA) (Frail et al. 1999a) since about 7.2 hours after 
the burst and referred as VLA J232221. 5+771543. The radio 
light curve shows a gradual rise to a plateau followed by a 
decline until below detectability after about 60 days. There 
are some large variations in these data, which is believed to 
be caused by interstellar scattering and scintillation (ISS; 
Frail et al. 1999a). 

As discussed by Dai & Lu (1999b), the steep decays 
of the X-ray and optical afterglows of GRB 980519 can be 
attributed to the shock evolution into the sub-relativistic 
phase 8 hours after the burst as the result of the dense cir- 
cumburst medium. During such a sub-relativistic expansion 
phase, the hydrodynamics of the shocked shell is described 
by the self-similar Sedov-von Neumann- Taylor solution. The 
shell radius and its velocity scale with time as r cx t 2 J° 

and f3 oc t^^ , where t® denotes the time measured in the 
observer frame. Then, we obtain the synchrotron peak fre- 
quency v m oc i^ 3 , the cooling frequency v c oc t^ 1 '' , the 
peak flux F Vm oc t 3 J° and the self-absorption frequency 

3p-2 

for the case of v a > ^ m (Dai & Lu 1999b). 



(1) 



{ iya/v m )-<*-W{ v / Va )*l*F Vm ocv^t 1 ^ 10 , 

if Um < V < V a \ 

{v/v m )-<*-WF Vm oc v-iP-V/^ 1 - 15 ^ 10 , 

if v a < v < v c ; 
{vc/vm)- {p - 1)/2 (v/v c )- p/2 F„ m oc u-^/H%-^' 2 , 
if v > v c . 



If the observed optical afterglow was emitted by slow-cooling 
electrons while the X-ray afterglow from fast-cooling electrons 
and if p S3 2.8, then according to Eq.(l), the decay index 
a R = (21 - 15p)/10 m -2.1 and a x = (4 - 3p)/2 »J -2.2, in 
excellent agreement with observations. Also, the model spectral 
index at the optical to X-ray band, (i = — (p— 1)/2 —0.9 is quite 
consistent with the observed one — 1.05±0.10. Furthermore, from 
the information of X-ray and optical afterglows, Dai & Lu (1999b) 
have inferred the physical parameters of this burst as follows: 

E~ 0.3 x 10 52 crg, e e ~ 0.16, 
n ~ 3 X 10 5 cm- 3 , z ~ 0.55, 



e B ~ 2.8 X 10~ 4 , 



(2) 



where E is shock energy, z is the redshift of the burst and e e and 
e b are the electron and magnetic energy fractions of the shocked 
medium, respectively. 

After the 60-day radio observational data being published, 
we promptly checked the dense medium model, and found that 
the asymptotic analysis can approximately describe the radio 
behavior. The analysis is as follows: adopting the inferred val- 
ues of the physical parameters in Eq.(2), the detected frequency 
v* = 8.46GHz equals to u a at about day 12; thus, according to 
Eq.(l), we expect that before this time the radio emissions rise as 
t^ 1 and then decay as tZ 2 ' 1 after the self-absorption frequency 
v a falls below v*. This simple asymptotic solution agrees quali- 
tatively with observations, as showed in Fig. 1 by the dotted line 
This preliminary analysis stimulated us to fit the radio data with 
a more detailed model by taking into account the trans-relativistic 
shock hydrodynamics and the strict self-absorption effects of the 
synchrotron radiation. 



3 TRANS-RELATIVISTIC SHOCK 

HYDRODYNAMICS, SELF-ABSORPTION 
EFFECT AND THE FITTING OF THE 
RADIO DATA 

We consider an instantaneous release of a large amount of energy 
E in a constant density external medium. The energy released 
drives in the medium a shock wave, whose dynamic evolution 
from the relativistic to sub-relativistic phase can be described 
approximately in the following way. 

Let r be the shock radius, 7 and Y be, respectively, the 
Lorentz factors of the shell and the shock front, f) be the velocity 
of the shock front. As usual, the shock expansion is assumed to 
be adiabatic, during which the energy is conserved, and we have 
(Blandford & McKee 1976) 



-TTaf3 2 r 2 r 3 nm p c 2 = E, 



(3) 



where u is a coefficient: a — » 0.35 when (i — > 1 and a — » 0.73 
when /3 — > 0. As Huang, Dai & Lu (1998), we use an approximate 
expression for a: a = 0.73 — 0.38/3. 

The radius of the shock wave evolves as (Huang, Dai & Lu 
1998) 



v a oc t„ 



i!l= / 3 C7 ( 7 + v / 7 2-l)/(l + z ). 



Now, the derived spectra and light curves are 



(4) 
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and the jump conditions of the shock are given by (Blandford & 
McKee 1976) 



/ 77 + 1 , 77 + 1 2 
-n, e = — : —(7 — \)nm v c , 



7-1 



7-1 



(7 + 1) [7(7- 1) + 1] 2 
7(2- 7 )(7-l)+2 



(5) 
(6) 



where e' and n' are the energy and the number densities of 
the shell in its comoving frame and 7 is the adiabatic index, 
which equals 4/3 for ultra-relativistic shocks and 5/3 for sub- 
rclativistic shocks. A simple interpolation between these two lim- 
its 7 = 4 ^ 1 gives a valid approximation for trans-relativistic 
shocks (Dai, Huang & Lu 1999). Using the above equations, we 
can now numerically obtain the evolution of r(%) and 7(t©) in 
the trans-relativistic stage, given proper initial conditions. 

As usual, we assume that the distribution of relativistic 
electrons with the Lorentz factor 7 e takes a power-law form 
with the number density given by n(-y e )d'y e = C-f e ~ p d'y e above 
a low limit 7 m i n , which is determined by the shock velocity: 
Is — 21H1p^ _ 1) Also, the energy densities of elec- 



7 



(p-l) 

trons and magnetic fields arc assumed to be proportional to the 
total energy density e' in the comving frame as U' e = e e e' and 
B' ± = (Swege') 1 ^ 2 . Thus, from the standard theory of the syn- 
chrotron radiation (Rybicki & Lightman 1979; Li & Chevalier 
1999), we have the expressions of the effective optical depth and 
the self-absorbed flux 

r„, = J^^(^^r p/2 F2(^)CB'^ + ^Ar', (7) 



8-nmv' 2 2mc 2 



Anmc 



3q 



C = (p-l)n' 7 ^, Ar' = r/ V , 



(8) 



F v = (1+ z )D 3 tt( — ) 



"(- 



3qB\ 



(l-e-),(9) 



where F\(x), F2{x) are defined by Eq.(5) in Li & Chevalier (1999), 
m and q denote the mass and charge of the electron, d^ is the 
luminosity distance of the burst, assuming a flat Friedman uni- 
verse with Ho = 65 km s^ 1 Mpc" 1 and r) ~ 10, characterizing 
the width of the shock shell. Here D = 1/(7(1 — 0)\ describes the 
relativistic effect, and v relates to the corresponding frequency v' 
in the comoving frame by v = Dv 1 /(l + z). 

Using the above full set of equations, we computed the ra- 
dio flux at the frequency = 8.46GHz and plotted the model 
fit in Fig. 1 as the solid line. We find that the following combi- 
nation of the parameters fits almost all valid data rather well: 
E ~ 0.8 X 10 52 erg, t e ~ 0.2, e B ~ 1 x 10~ 4 , n~lx 10 5 cm" 3 , 
z ~ 0.55 and the Lorentz factor at the initial time % = 1/3 days, 
7 = 1.2 (/3 ~ 0.55). We stress that these parameters are in ex- 
cellent agreement with those inferred independently from the X- 
ray and optical afterglows (Dai & Lu 1999b), as listed in Eq.(2). 
Clearly, there are some large amplitude variations in the observed 
light curve (e.g. about 18 days after the burst), which is believed 
to be caused by diffractivc scintillation. Also plotted in Fig. 1 is 
the fit (dashed line) computed with the sub-relativistic model as 
presented in the appendix of Frail, Waxman & Kulkarni (1999). 
The fit with this model was obtained by adopting the initial con- 
ditions as: to = 1/3 days and 7*0 = 1.6 X 10 16 cm and the parameter 
values (e e , eg, E, z and n ) the same as those in Dai & Lu (1999b). 
Comparing the trans-relativistic model with sub-rclativistic one, 
we can easily see that the relativistic effect (as characterized by 
D = l/7(l — /3) in Eq.(9)) flattens the rising phase at earlier 
time, making the trans-relativistic model agree better with the 
observations, while at the later time both the fitting curves trend 
towards the asymptotic solution (i.e. F v oc tZ 2 ' 1 ). 



3. 
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Figure 1. Model fits of the radio light curve at 8.46 GHz. De- 
tections and upper limits for the non-detections, taken from Frail 
et al. (1999a), are indicated by the filled squares and arrows, re- 
spectively. The dotted line represents the expected asymptotic 
solution of the radio behavior according to the dense medium 
model in Dai & Lu (1999b). The solid and dashed lines repre- 
sent the fits with the tran-relativistic and sub-relativistic models, 
respectively. See text for additional details. 



4 DUST SUBLIMATION AND OPTICAL 
EXTINCTION BY DENSE MEDIUM 

One may ask whether the dense circumburst medium may cause 
large extinction in the optical afterglow of GRB 980519. A crude 
estimate is as follows. At the time that the blast wave transits to 
the sub-relativistic stage (% ~ ^days, (3 ~ 0.55, 7 ~ 1.2), from 
Eq.(3) we derived its radius to be r ~ 2.1 X 10 16 cm. Therefore the 
characteristic column density through the medium into which the 
blast wave is expanding is about nr ~ 2.1 X 10 21 cm -2 with the 
corresponding Ay of 1.3 magnitudes in the rest frame of the ab- 
sorber. This column density is comparable to (but slightly larger 
than) the Galactic 21cm column density (~ 1.74 X 10 21 cm~ 2 , 
Halpern ct al. 1999) in the direction of the GRB 980519. With 
Ax scales linearly with i , the absorption in our observed B band 
for this optical transient is about 1.3(l+z) ~ 2 magnitudes, twice 
the value adopted by Halpern et al. (1999) for correction of the 
relative extinction. 

In the above estimate, we have made a questionable assump- 
tion, i.e. the dense medium around the burst has a standard gas- 
to-dust ratio. However, this may be not realistic, considering that 
the dust around the burst can be destroyed due to sublimation out 
to an appreciable distance (~ a few pc) by the prompt optical-UV 
flash (Waxman & Draine 1999; hereafter WG99), accompanying 
the prompt burst. Below we will give an estimation of the de- 
struction radius for dense medium case, following WG99. 

The prompt optical flash detected accompany- 
ing GRB 990123 (Akcrlof ct al. 1999) suggests that, at least for 
some GRBs, 7-ray emission is accompanied by prompt optical- 
UV radiation with luminosity in the 1-7. 5eV range of the order 
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of 10 49 (^-)erg/s for typical GRB parameters, where AS7 is the 
solid angle into which 7-ray and optical-UV emission is beamed 
(WG99). The most natural explanation of this flash is emission 
from a reverse shock propagating into the fireball ejecta shortly 
after it interacts with the surrounding gas (Sari & Piran 1999; 
Meszaros & Rees 1999). As for GRB 980519 with the parameter 
values as e e ~ 0.2, e B ~ 10~ 4 , n ~ 10 5 cm~ 3 , E ~ 8 X 10 51 erg 
and the burst duration At ~ 70s, we derive the luminosity in the 
1-7.5 eV range to be about L1-7.5 ~ 5 X 10 48 erg/s (Here we have 
assumed that electron and magnetic field energy fractions in the 
reverse shock are similar to those in the forward shock; Wang, 
Dai & Lu 1999c). The condition for the grain to be completely 
sublimed during the prompt flash time is 

T >T C ~ 2300K[1 + 0.0331n( °~ 5 )], (10) 

At/ 10s 

where T is the grain temperature, determined by Eq.(8) of WG99, 
and a = a— 5 X 10~ 5 cm is the radius of the dust grain. Then, 
according to Eq.(17) of WG99, the radius out to which the prompt 
flash can heat grains to the critical temperature T c is 

R c ~ 3.7xlQ 19 ( °^ L49(1 + ala - 5) ) 1 / 2 cm ~ 2.7x 10 19 cm,(ll) 

where Qjjy is the absorption efficiency factor of the optical-UV 
flash and can be assumed to be near one for grain radii a > 
10~ 5 expected in dense medium. Thus, we can safely say that the 
extinction due to the circumburst dense medium is not important 
if the size of the dense medium is shorter than R c and we think 
that this condition is reasonable for GRB 980519, and also for 
GRB 990123 (Dai & Lu 1999a). 



5 DISCUSSIONS AND CONCLUSIONS 

The detection of strong diffractivc scintillation requires that the 
angular size of the source VLA J232221. 5+771543 should be 
less than l^tarcsec even 15 days after the burst (Frail et al. 
1999a), otherwise the fluctuations would be suppressed. This 
small inferred size is not consistent with the spherical, homo- 
geneous model with a normal density (n ~ lem -3 ), but can 
marginally be consistent with the jet model and the wind- 
shaped circumburst medium model. We note that the dense 
medium model may have a potential advantage for this require- 
ment, simply because that the shock will be quickly deceler- 
ated to a sub-relativistic velocity and therefore have a shorter 
shock radius. This can be clearly seen from the following com- 
parison: s re i oc J* o ® 27 2 cdt(l/-y(t®)) for the relativistic case 

while 9 a ^ auo oc j*® 0cdt(f3 < 1) for the sub-rclativistic one. 

Here, for the afterglow of GRB 980519, we assume that before 

tg5 ~ l/3days, the shock is adiabatic and relativistic; thus the 

shock radius is r(t) ~ (17 'Et® / 'iimi p nc{l + z)) 1/4 ~ 0.78 X 
10 16 (£/3x 10 51 crg) 1 /4( n / 3x i 5 cm -3)-i/4( tffl /i days )i/4 [( - 1 + 

z)/1.55] -1 / 4 cm (Sari, Piran & Narayan 1998). Then, the shock 

radius should follow the Sedov-von Neumann- Taylor self-similar 

2/5 

solution as r(t®) oc tJ ■ Thus, we obtained the angular size of the 
afterglow 9 3 ~ O.8/iarcsec(t0/15days) 2 / 5 . The agreement would 
be even improved if we note that at the beginning of the assumed 
sub-relativistic stage, the radius should increase more slowly than 
the self-similar solution in the trans-relativistic regime. 

The strong modulations caused by scintillation also make 
the estimate of the spectral slope in the radio band less accu- 
rate. The averaging value of the spectral slope from day 12 on 
is (3 ~ —0.45 ±0.6 (where F v oc u^), implying that the time 
averaged self-absorption frequency v a is between 1.43 GHz and 



4.86 GHz (Frai et al. 1999b). In our model, the time (t e ~ 
12days) when the fitting curve begins to decline corresponds to 
v a = v* = 8.46GHz. Since the self-absorption frequency decays 
as v a oc ( e -(3p-2)/(p+4) ^ t -0-94 for p = 2 8] we expect that 

u a shifted quickly below 4.86 GHz at day 21, but was above 1.43 
GHz over all the detecting time, which is in reasonable agreement 
with the observations. 

The radio afterglow of GRB 970508, the longest light curve 
(450-day) obtained by far, exhibited different behavior from 
GRB 980519 (Frail et al. 1997; Waxman, Kulkarni & Frail 1998; 
Frail, Waxman & Kulkarni 1999). From the spectral and temporal 
radio behavoir, Frail, Waxman & Kulkarni (1999) inferred that 
the fireball has undergone a transition to sub-rclativistic expan- 
sion at t ~ 100 days, consistent with the inferred low ambient den- 
sity n ~ lcm~ 3 (but also see Chevalier & Li 1999a). On the other 
hand, some radio afterglows (e.g. GRB 990510, GRB 981226; Frail 
et al. 1999b) show similar behaviors to GRB 980519, that is, they 
exhibit a slow rise to the maximum for a relatively short time and 
then a fast decline until below detectability. It is likely that the 
shocks of these bursts entered into the sub-relativistic stage after 
a short relativistic one and our above model can also describe 
their radio afterglows. Harrison et al. (1999) had interpreted the 
broad-band lightcurve break in the afterglows of GRB 990510 as 
due to a jet-like outflow. We speculate that another possible ex- 
planation is that the shock had entered into the sub-relativistic 
stage after ~ ldays as the result of the combination of the dense 
medium and jet effects (Wang et al. 1999b), the latter of which 
may be real in consideration of the large inferred isotropic energy. 
The radio afterglow of GRB 990123 is unique for its "flare" be- 
havior (Kulkarni et al. 1999b), whose most natural explanation is 
that it arises from the reverse shock, as evidenced by the prompt 
optical flash (Sari & Piran 1999). Our preliminary computation 
(using the trans-relativistic model) shows that the radio emission 
from the forward shock in the dense medium model is significantly 
lower than that from the reverse shock and declines quickly after 
the peak time, if a jet-like outflow with an opening angle ~ 0.2, 
as required by the "energy crisis" of this burst, is invoked (Wang 
et al. 1999b). Moreover, the fast decline of the radio emission 
from the forward shock, which is caused by the deceleration of 
the shock in the sub-relativistic stage, can be consistent with the 
non-detection even 3 days after the burst. 

In summary, we argue that the dense medium model, which 
has interpreted the optical to X-ray afterglows of GRB 980519 
quite well, can also account for the radio afterglow excellently. The 
circumburst environment can affect the evolution of GRBs after- 
glows significantly (Meszaros, Rees & Wijers 1998; Panaitescu, 
Meszaros & Rees 1998; Wang et al. 1999a). For a low (n ~ 
lem -3 ), homogeneous density environment, the shock waves stay 
at the relativistic shock stage for quite a long time, while for the 
dense medium case, the shock wave quickly enters into the sub- 
relativistic stage. Recently, a generic dynamic model for the evolu- 
tion of shocks from ultra-relativistic phase to the sub-relativistic 
one has been also developed by Huang et al. (1999a). The af- 
terglows of the optically thin radiation (e.g. optical and X-rays) 
from the shock at the sub-relativistic stage decays more rapidly 
than at the relativistic one. As for the radio afterglow (usually 
% < v a at the sub-rclativistic stage for this model), the dense 
medium model predicts a slow rise (y* < v a ), followed by a round 
peak and a late steep decline (v t > u a ), trending towards the be- 
havior of the optical and X-ray afterglows. Clearly, this behavior 
is different from the jet model in the early epoch. But it is some- 
what similar to the wind model, making it difficult to distinguish 
between them through the radio observations. 
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